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ABSTRACT 


An  approximate  theory  for  cavitation  in  pump  impellers  is 
developed  on  the  basis  of  a  one-dimensional  potential  flow 
model.   The  theory  incorporates  the  unsteady  term  occurring 
in  Bernoulli's  equation  and  accounts  for  blockage  effects  due 
to  vane  thickness.   A  comparison  of  the  theoretical  results 
with  experimental  data  from  mixed  flow  impeller  tests  shows 
good  agreement  near  the  condition  of  zero  incidence.   On  the 
basis  of  the  theoretical  results  a  mechanism  of  cavitation 
breakdown  is  proposed. 


NOMENCLATURE 


2 
A  =  pipe  area  upstream  of  impeller,  ft. 

2 
A.  =  area  of  flow  passage  at  the  impeller  eye  without  vanes,  ft. 

B  =  vane  blockage  coefficient,  dimensionless 
v 

C  =  -; — - —  =  vane  coefficient,  dimensionless 
v    1-B 
v 

C  =  shroud  coefficient,  dimensionless 
s 

C  =  meridional  velocity,  ft/sec. 
m 

D  =  impeller  diameter,  ft. 

g  =  gravitational  constant,  ft/sec. 

H  =  net  positive  suction  head,  ft. 
sv      r 

K  =  cavitation  index,  dimensionless 

N  =  impeller  rotative  speed,  rpm 

2 
P  =  pressure,  lb/ft. 

2 


P  =  vapor  pressure,  lb/ft 
v 

q  =  velocity  vector,  ft/sec. 

3 
Q  =  flow  rate,  ft  /sec. 

t  =  time,  seconds 

u  =  impeller  wheel  speed,  ft/sec. 

Z  =  number  of  impeller  vanes 

a  =  relative  flow  angle,  radians 

8  =  vane  angle,  radians 

.  .  -1 

5  =  vorticity  vector,  sec. 

x  =  blade  thickness,  in. 

<f>   =  C  /u,  flow  coefficient,  dimensionless 
m 

2 
$  =  velocity  potential,  ft  /sec. 

SUBSCRIPTS 

0  =  reference  station  upstream  of  the  impeller 

1  =  reference  station  outside  vane  row  at  impeller  inlet 

2  =  reference  station  inside  vane  row  at  impeller  inlet 


A  THEORY  FOR  CAVITATION  IN  IMPELLERS 

INTRODUCTION 

It  has  long  been  recognized  that  the  flow  in  turbomachinery 
is  by  nature  unsteady  when  viewed  from  a  fixed  reference  frame. 
Dean  [l]  pointed  out  the  necessity  of  unsteadiness  by  showing 
that  Bernoulli's  equation,  applied  between  the  inlet  and  out- 
let of  a  pump,  denies  the  addition  of  energy  to  the  fluid 
unless  the  unsteady  terms  are  included.   In  other  words,  energy 
cannot  be  added  to  a  steady,  irrotational  flow  of  a  friction- 
less  fluid;  only  by  making  the  boundary  conditions  unsteady  is 
it  possible  to  add  energy  to  the  fluid. 

As  a  consequence  of  the  inherent  unsteady  flow  occurring  in 
impellers  when  viewed  from  a  fixed  reference  frame,  it  is 
necessary  to  account  for  the  unsteady  effects  in  any  theoret- 
ical attempt  to  predict  cavitation  or  minimum  pressure  coeffi- 
cient.  Accordingly,  the  object  of  this  paper  is  to  develop  a 
theory  for  cavitation  in  pump  impellers  which  accounts  for  such 
effects.   Bernoulli's  equation  is  applied  to  the  flow  in  the 
inlet  region  of  an  impeller  and  the  unsteady  term  occurring 
therein  is  retained.   An  analysis  is  then  developed  to  approxi- 
mate this  term  so  that  the  final  result  yields  the  cavitation 
index  as  a  function  of  the  flow  coefficient. 

One  of  the  most  surprising  results  of  this  analysis  is  that 
the  unsteady  term  in  Bernoulli's  equation  is  not  of  secondary 
importance  but  has  a  major  influence  on  the  pressure  in  the 
inlet  region  at  values  of  the  flow  coefficient  on  either  side 
of  the  zero  incidence  condition.   Moreover, on  the  basis  of  the 


theoretical  results  obtained,  a  mechanism  for  cavitation 
breakdown  at  positive  and  negative  incidence  conditions  is 
proposed . 

The  underpressure  in  the  eye  of  an  impeller  has  been 
considered  to  be  dependent  on  two  factors,  the  effect  of 
the  shroud  and  the  impeller  vanes.   To  obtain  the  total 
maximum  underpressure  in  the  eye,  the  vane  pressure  drop 
is  added  to  the  pressure  drop  caused  by  the  shroud.   Reck- 
oning in  this  manner,  Spannhake  and  Wislicenus  developed  an 
expression  for  the  cavitation  index  in  the  form 
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where  Ki  and  K2  are  shroud  and  vane  coefficients,  respectively, 
Q  denotes  the  flow  rate  and  D  denotes  the  impeller  diameter. 
Oongwer    [-J  was  a°le  to  draw  a  mean  line  through  his  experi- 
mental data  by  properly  selecting  values  for  K}  and  K2  and, 
in  terms  of  the  notation  of  this  paper,  his  resulting  empirical 
expression  for  cavitation  index  at  head  breakoff  versus  flow 
coefficient  at  zero  incidence  is 


V1-4  +  -085  1+  -v 


Although  this  empirical  equation  represents  Gongwer's  experi- 
mental data  for  impellers  having  similar  eye  designs  fairly  well 
it  gives  little  insight  into  the  mechanism  involved  and  does  not 
reflect  the  effect  of  such  features  as  vane  thickness  or  number 
and  non-zero  incidence  angle  of  the  flow. 


An  analytical  method  to  calculate  the  complete  flow  field 
in  an  impeller  is  outlined  in  Ref.  [7]  .   This  stream  filament 
technique  was  programmed  on  a  large  digital  computer  and  used 
to  design  a  number  of  mixed-flow  impellers  described  in  Ref. 
|_6J.   A  comparison  of  the  test  results  of  these  impellers 
with  the  theory  showed  the  theory  to  have  the  proper  general 
trends  and  give  reasonably  accurate  results  in  the  portion  of 
the  impeller  where  channel  flow  exists.   However,  in  the  inlet 
region  of  the  impeller  the  theory  gave  pressures  which  were 
consistently  too  high,  making  it  of  little  value  for  the  pre- 
diction of  cavitation. 

Cavitation  in  impellers  and  inducers  has  generally  been  con- 
sidered to  be  described  by  separated  cavity  flow  at  the  leading 
edge  of  the  vanes,  similar  to  that  which  is  observed  in  steady 
flow  past  a  single  foil.   Analytical  methods  to  describe  this 
kind  of  cavitation  breakdown  [.2j  ,  [3]  ,  [4j  have  been  based  on 
the  model  of  a  free-streamline  flow  through  a  cascade  of  flat 
plate  hydrofoils.   While  this  model  is  no  doubt  a  reasonable 
representation  of  cavitating  flow  at  large  positive  angles  of 
incidence,  it  is  of  limited  value  near  zero  and  at  negative 
angles  of  incidence  (overflows).   The  present  theoretical 
findings  indicate  that  at  negative  incidence  a  favorable  pres- 
sure gradient  exists  and  there  is  no  tendency  for  a  separated 
region  to  exist. 

One  of  the  first  rather  comprehensive  contributions  to  the 
literature  on  pump  cavitation  was  published  by  Gongwer  [2] . 
A  series  of  impellers  of  approximately  the  same  eye  design 
with  widely  differing  outside  diameters  were  tested  with  the 


result  that  the  effect  of  the  differing  outside  diameter  had 
no  effect  on  cavitation.   Thus,  Gongwer  concluded  that  the 
cavitation  performance  of  an  impeller  is  only  dependent  on 
the  eye  design  and  the  angle  of  incidence  that  the  fluid  makes 
with  the  impeller  vanes. 

A  simplified  flow  model  advanced  by  Stanitz  [5]  to  describe 
the  effect  of  vane  blockage  at  zero  incidence  was  applied  by 
Wood,  Murphy  and  Farquhar  [pi     to  their  experimental  results 
from  mixed-flow  pump  tests.   This  model  accounts  for  the 
reduced  area  caused  by  blockage  effects  of  the  vanes  according 
to  the  one-dimensional  continuity  equation  and  relates  the 
pressure  at  a  station  upstream  of  the  impeller  to  the  pressure 
at  a  station  just  inside  the  impeller  according  to  the  steady- 
flow  form  of  Bernoulli's  equation.   The  theoretical  approach 
presented  herein  is  essentially  an  extension  of  Stanitz's 
model  to  include  the  time  derivative  term  in  Bernoulli's 
equation . 

The  results  of  two  outstanding  experimental  investigations 
of  cavitation  in  mixed  flow  impellers  are  presented  in  the 
companion  papers,  Ref.  [6j  and  [8j .   The  tests  carried  out 
under  carefully  controlled  conditions  showed  the  effect  of 
vane  incidence  angle,  thickness  and  number  as  well  as  speed 
effects.   Since  the  geometrical  characteristics  of  the  series 
of  impellers  tested  are  well-recorded .  the  results  were  found 
to  be  very  valuable  for  comparison  with  the  theoretical  result 
developed  in  this  report. 


THEORY 

Two  probable  types  of  cavitation  are  hypothesized:  (a) 
localized  cavitation  at  the  leading  edge  of  the  vanes,  (b) 
gross  cavitation  behind  the  leading  edge,  which  will  be 
considered  to  be  the  true  cavi ta t ion, giving  rise  to  a 
significant  reduction  in  performance.   While  the  blunt  lead- 
ing edge  part  of  the  vanes  is  no  doubt  the  critical  point  in 
the  impeller,  it  is  believed  that  cavitation  here  would  be 
only  a  local  phenomenon  and  probably  would  not  result  in  a 
reduction  in  performance.   Thus,  we  will  consider  the  point 
just  behind  the  leading  edge  as  the  critical  point  for  gross 
cavitation  and  attempt  a  theoretical  prediction  of  the  pres- 
sure or  cavitation  index  at  this  point. 

In  essence,  the  basic  problem  is  to  calculate  the  minimum 
pressure  coefficient  on  the  impeller  under  unsteady  flow 
conditions.   Assuming  the  flow  to  be  incompressible  and 
inviscid,  Bernoulli's  equation  applied  between  some  point 
in  the  inlet  conduit, far  enough  removed  from  the  impeller 
eye  so  that  the  flow  can  be  considered  locally  steady,  and 
some  point  within  the  impeller,  gives: 
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In  this  equation  q  is  the  local  velocity  measured  with  respect 
to  a  fixed  reference  frame.   The  subscript  m  denotes  "meridonal' 
and  "0"  corresponds  to  the  upstream  reference  point.  The 
velocity  potential  $  is  defined  as  q  =  g  r  a  d  $  =  V  $  .  Neglecting 
the  elevation  term  and  adding  and  subtracting  the  fluid  vapor 


pressure,  P  ,  gives 
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The  classical  theory  of  cavitation  states  that  cavitation 
occurs  when  the  local  minimum  pressure  just  equals  the  vapor 
pressure  of  the  liquid.   It  is  well-known,  however,  that  in 
practice  this  simple  concept  is  not  always  valid.   Departures 
of  the  cavitation  index  from  values  calculated  on  the  basis 
of  this  hypothesis  are  caused  by  so-called  scale  effects 
which  arise  from  changes  in  size,  velocity  and  fluid  properties 
Fortunately,  however,  experimental  evidence  indicates  [9]  that 
in  the  absence  of  thermodynamic  effects  the  scale  effects 
decrease  with  size  and  velocity  so  that  for  the  prototype  ,  this 
hypothesis  yields  valid  results.   Thus,  proceeding  on  the  basis 
that  the  size  and  speed  of  the  impeller  are  sufficiently  large 
so  that  the  condition  of  cavitation  is  reached  when  P=P  , 


equation  (1)  may  be  written 
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where 
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mo 


is  defined  as  the  cavitation  index.   The  cavitation  index  in 
this  formulation  is  exactly  equal  to  one  minus  the  minimum 
pressure  coefficient.   Unfortunately,  this  definition  differs 
by  a  factor  of  1.0  with  the  index,  a,  which  is  in  common  use 
for  uniform  flow  past  bodies  such  as  head  forms  and  hydrofoils. 


It  is,  however,  consistent  with  the  definition  generally  used 
for  turbomachinery . 

According  to  equation  (2)  it  is  evident  that  in  order  to 

calculate  the  cavitation  index,  it  is  necessary  to  evaluate 


both  the  absolute  local  velocity,  q,and 


3  $ 
3t 


at  the  critical 


location  on  the  impeller.   As  a  first  approximation,  q  may  be 
taken  to  be  the  meridional  velocity  at  the  point  in  question, 
which  will  be  dependent  on  the  shroud  and  hub  configuration, 
corrected  for  vane  blockage.   The  effect  of  vane  thickness 
can  be  approximated  by  accounting  for  the  reduction  in  flow 
area  presented  by  the  vanes  according  to  one-dimensional  flow 
pr inc  ip les . 

The  term  3$/3t  can  also  be  approximated  once  q  is  known 
provided  the  kind  of  unsteady  flow  involved  is  recognized. 
That  is,  although  the  flow  is  unsteady  when  viewed  from  a 
fixed  reference  frame,  it  is  steady  in  a  reference  frame  rota- 
ting with  the  impeller.   As  a  consequence,  we  may  imagine  a 
velocity  potential  field  fixed  to  and  rotating  with  the  impeller 
The  term  3$/3t,  which  denotes  the  time  derivative  of  $  measured 
in  the  fixed  reference  frame,  may  be  calculated  by  observing 
the  change  in  $  as  the  $-field  rotates  past  a  point  fixed  in 
space  during  a  small  time  interval. 

If  the  vane  thickness  was  zero  and  the  vanes  were  aligned 
with  the  flow,  the  surfaces  of  constant  $  would  represent  disks 
or  surfaces  axisymmetric  with  respect  to  the  impeller  axis. 


In  this  case 


3  $ 
3t 


would  vanish  since  an  observer  fixed  in  space 


would  observe  only  one  value  of  0  as  these  disks  are  rotated. 
However,  when  the  vanes  have  a  finite  thickness  or  are  operated 


at  some  angle  of  incidence,  the  flow  is  disturbed  and  the 
surfaces  of  constant  $  are  distorted  so  that  they  are  no 


longer  axisymmetr ic .   In  this  case 


9$ 


3t 


would  no  longer  be 


zero  . 

Figure  1  shows  a  schematic  of  the  unperturbed  flow  where 

the  vane  thickness  is  zero  and  the  vanes  are  perfectly  aligned 

with  the  flow.   In  this  case  C   =  C  0  and  the  velocity  potential 

ml   m2 

lines  are  aligned  with  the  vane  velocity  vector  as  indicated. 
For  the  flow  depicted  here,  3$/8t=0  because  the  lines  corres- 
ponding to  constant  values  of  $  move  in  the  direction  of  their 
orientation.   An  observer  fixed  in  space  would  observe  the  same 
value  of  $  as  time  progresses  and  the  velocity  potential  field 
rotates.   For  this  special  unperturbed  condition  the  cavitation 
index  as  given  by  equation  (2)  reduces  to 

K  =(C  ,  ,C   )2 

o    ml /  mo 


This  relationship  which  represents  the  lower  limit  of  K  reflects 

r  r  o 

only  the  effect  of  the  shroud  on  cavitation  index  since  the  vane 

have  no  thickness  and  are  aligned  with  the  flow. 

Figure  2  illustrates  the  case  of  finite  vane  thickness,  i.e., 

a  small  perturbation  of  the  flow  depicted  in  Figure  1.   For  this 

case  it  is  evident  that  C   >C  ,  where  C  „  may  be  related  to  C  . 

m2   ml  m2  ml 

approximately  by  the  one-dimensional  continuity  relationship 
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where  the  vane  blockage  coefficient,  B  ,  represents  the  ratio 
of  the  vane  frontal  area  to  the  area  of  the  passage  without 


vanes.   For  purposes  of  evaluating  B   the  projected  vane 
frontal  width  is  taken  as  x/sinB  where  t  is  the  vane  thick- 
ness and  6  is  the  local  vane  angle  indicated  in  Figure  2. 

3  $ 
Since  C  „  is  greater  than  C  .  and  C  =  - — ,  the  velocity 
mz  ml       m   3x  J 

potential  lines  must  be  more  closely  spaced  in  the  axial 

direction  for  finite  vane  thickness.   Nevertheless,  the 

basic  kinematic  boundary  condition  on  the  surface  of  the  vane, 

given  by 

a  <*    -  - 

(4) 
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where  n  denotes  the  unit  normal  vector  to  the  blade  and  u 
denotes  the  local  vane  velocity  vector,  must  still  be  satisfied 
as  it  was  for  the  case  of  zero  vane  thickness.   That  is,  in 
either  case  the  normal  derivative  of  the  velocity  potential 
must  satisfy  the  condition: 


3<3> 
9n 
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From  this  condition,  then,  the  new   1  and   2  lines  for  the 

perturbated,  flow  having  the  same  values  as  those  indicated  in 

Figure  1, can  be  drawn.   Thus,  in  Figure  2  we  have  represented 

by  dashed  lines  the  new   1  and   2  lines  which  are  more  closely 

spaced  in  the  x-direction  as  required  since  C   >C  .  but  still 

m2   ml 

3  $ 
having  the  same  value  of  7 —  on  the  vane  surface.   That  is,  to 

satisfy  this  normal  derivative  requirement  the  new.$  lines  are 

constructed  so  as  to  intersect  points  B  and  C  in  Figure  2. 

n,  3  $  3  $     , 

To  evaluate  - — ,   A$  can  be  written  as  A$=3>  -<P  =  - —  Ax   where 
3  t  2   1   dx 

Ax*  is  the  axial  distance  indicated  in  Figure  2.   The  corres- 

Ac 

ponding  At  is  At=  where  Ac  is  the  distance  indicated  in 
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Figure  2  and  u  is  the  vane  velocity.   Thus, 


21    A$    3x 
3t    At 


Ax' 


3$   Ax' 


Ac         3x   Ac 


(5) 


but,  as  indicated  in  Figure  2, 
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3  $ 
for  small  values  of  0.  Thus,  from  equation  (5)  and  (6)  -r—   can 

a  t 

be    written , 
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3  $ 
where  - —  has  been  replaced  by  C  „. 
3x  m2 
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By  definition  of  the  velocity  potential,  - —  is  equal  to  the 

3n 

velocity  of  the  fluid  in  a  direction  normal  to  the  vane  surface 

From  Figure  2  it  is  evident  that  this  can  be  written  as 
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Considering  the  geometry  of  the  triangle  (BCD),  indicated  in 
Figure  2f   it  is  evident  that 


a=  Ax  tan 


(9) 


and 


Ax-Ax'=a  tan0  -  a0  (10) 

Eliminating  a  between  equations  (9)  and  (10)  and  substituting 
this  result  into  equation  (8)  gives 

(11) 


- —  =  -  t —  (1-OtanB)  cos 
3n      3x 


3  $ 
Substituting  this  expression  for  - —  into  the  kinematic  boundary 

d  n 

condition  given  in  equation  (4)  yields 


0  =cot8- 


C  2 

m 


where  ,  again  , 


3$ 

3x 


has  been  replaced  by  C  n.       Substituting  this 


expression  into  equation  (7)  yields  the  following  expression 


for 


3$ 


~—   in  terms  of  the  vane  velocity,  angle  and  meridional 

velocity : 

3$ 


d  *  /         U    \ 

-  -  u  cm2(cote-  — ) 

m2 


(12) 


Finally,  substituting  the  relationship  given  in  equation  (12) 

3  $ 
for  — —  and  using  the  approximation  q  =  C  0  in  equation  (2) 
d  t  m  Z 

yields  the  following  relationship  for  the  cavitation  index: 


K  = 
o 


mz 


mo 


+  2 


u m2 

C     C 

mo    m  1 


(cote-  =_) 

m2 


(13) 


Since,  in  general,  the  upstream  reference  velocity,  C   , 

mo 

and  the  local  shroud  velocity  for  zero  vane  thickness,  C  , , 

m  1 

differ  due  to  the  effect  of  the  shroud,  it  is  convenient  to 

introduce  C  ,  into  equation  (13)  to  obtain 
m  I 

K  JjZLJ   |^2_\  +  2  -4-   ^McotB-  ~±      Cm°  CmlA( 


o  \C    / \C  , 

\  mo  '  \  ra 1 


*    C 
o    ml 


C  .  C  0 
o   ml   m2 


where  d>  =C   /u  is  the  flow  coefficient  based  on  the  upstream 
o   mo 

meridional  velocity.   C  ,/C         can  be  called  a  shroud  coeffi- 

ml   mo 

cient  since,  from  a  one-dimensional  flow  point  of  view,  it 

represents  the  ratio  of  the  inlet  pipe  reference  area  to  the 

local  passage  area  if  vanes  are  not  present.  C  _/C  ,  can  be 
°  m2   ml 

called  a  vane  coefficient  since  it  represents  the  ratio  of  the 
local  passage  area  without  and  with  vanes,  respectively.  From 
one-dimensional  continuity  considerations,  as  given  in 
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equation  (3 )  , 


C  =  7 
v   C 


m2 


ml 


1-B 


=vane  coefficient 


(15) 


v 


and 


C  = 
s 


'ml 


=  shroud  efficient 


(16) 


mo 


where  B   is  the  vane  blockage  coefficient  and  A  /A.  is  the  area 
v  o   1 

of  the  reference  section   in  the  inlet  pipe  to  the  local  area  of 
the  passage  excluding  the  effect  of  vane  blockage.   Using  these 
definitions,  the  final  expression  for  the  cavitation  index 

becomes 

„       ~    2-  2  .  -  „      ..1      2   1 

(17) 


K  =C   C    +  2  C    cotB— r 
o   s   v        v        <p 


C   <f>  z 
s   o 


Equation  (17)  relates  the  cavitation  index  to  vane  and  shroud 
coefficients,  the  vane  angle  and  the  flow  coefficient.   The 
first  term  represents  the  dynamic  depression  associated  with  the 
velocity-squared  term  in  the  Bernoulli's  equation  while  the 
second  two  terms  together  represent  the  contribution  from  the 
unsteady  term.   The  equation  is  valid  at  any  radial  location  on 
the  vane  if  3  and<}>0  are  interpreted  as  the  local  vane  angle  and 
local  flow  coefficient,  respectively.   However,  numerical  evalua- 
tion of  the  equation  shows  that  the  maximum  value  of  K   occurs 

o 

at  the  vane  tip  as  might  be  expected.   Accordingly,  3  and  d> 

o 

should  be  interpreted  as  the  vane  setting  angle  and  ratio  of 
cmo  to  the  vane  wheel  speed  at  the  vane  tip,  respectively.   It 
may  be  noted  also  that,  although  the  "suction"  side  of  the  vane 
was  used  in  the  analysis,  the  resulting  equation  is  the  same 
regardless  of  which  side  is  employed.   It  is  only  necessary  that 


8  be  interpreted  as  the  angle  of  the  vane  face  in  question 
if  they  are  different. 

As  indicated  previously  the  cavitation  index  is,  by 
definition,  related  to  the  minimum  value  of  the  pressure 
coefficient  according  to 


K  - 1  —  C  / _  •  _ \ 

o     p  (.mm; 


where  C  ,  .  Nis  based  on  the  pressure  and  velocity  at  the  up- 
p  (mm )  ^  j  v 

stream  sub  "o"  reference  station.   A  value  of  K  =1.0  therefore 

o 

corresponds  to  the  condition  where  the  local  pressure  is  equal 
to  the  upstream  reference  pressure.   This,  of  course,  would 
correspond  to  the  minimum  possible  value  of  the  cavitation 
index  excluding  effects  of  both  the  vanes  and  shroud.   Values 

of  K   less  than  unity  would  simply  indicate  that  the  local 

o  '  v    J 

pressure  was  greater  than  the  reference  pressure  and,  these 

values  are  of  no  significants. 

The  condition  where  the  flow  just  inside  the  vanes  at  station 

"2"  is  aligned  with  the  vanes  is  of  particular  interest.   This 

condition  is  defined  by  <J>  =tanB  or,  equiva  len  t  ly ,  <$>    =tanB/C  C  . 

2  o         s  v 

Substituting  this  value  for  cj)   into  equation  (17)  yields  simply 


K  =C  2C  2 
o    S    V 


(18) 


because  the  second  and  third  terms,  corresponding  to  the  contri- 
bution of   3$/3t,  cancel  each  other.   This  result  is  essentially 

that  of  the  blockage  model  proposed  by  Stanitz  L5j  . 

Values  of  <t>       which  are  less  than  tan8/C  C   correspond  to 

0  S  V 

underflows  or  positive  incidence  at  the  vanes  and  the  contribu- 
tion from  3  $  /  3  t  in  equation  (17)  is  negative,  tending  to  decrease 
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K  .   This  effect,  interpreted  in  terms  of  local  pressure, 
o 

corresponds  to  an  increase  in  pressure  and  pressure  coef- 
ficient .  Although  the  measurements  [6J  show  that  the  average 
pressure  inside  the  vanes  does,in  fact, increase  for  the 
positive  incidence  condition,  it  is  also  clear  that  the 
measured  cavitation  index  increases,  indicating  a  local 
depression.   It  is  understandable  then  that  Equation  (17), 
which  is  based  on  one-dimensional  flow  considerations, 
disregards  the  local  anomolies  due  to  leading  edge  flow 
and  predicts  the  increase  in  pressure  at  stations  behind 
the  leading  edge.   This  positive  incidence  flow  is  dis- 
cussed in  more  detail  and  in  relation  to  experimental 
findings  in  a  later  section. 

Flow  at  values  of  d>   which  are  greater  than  tan3/C  C 

o  s  v 

corresponds  to  underflow  or  a  negative  incidence  at  the 
leading  edge  of  the  vanes.   For  this  condition  the  last 
two  terms  in  equation  (17),  which  represent  the  contribu- 
tion from  3$/3t,  give  a  positive  contribution  to  K  .   This 
effect,  when  interpreted  in  terms  of  pressure,  corresponds 
to  a  depression  and,  consequently,  reduced  cavitation  perfor- 
mance . 
COMPARISON  WITH  EXPERIMENTAL  RESULTS 

Figures  3,4,5  and  6  show  the  experimental  results  from 
mixed-flow  pump  tests  of  Wood,  Murphy  and  Farquhar  [6 J  and 
Wood  [8~\      along  with  the  present  theoretical  results  obtained 
from  equation  (17).   The  experimental  points  represent  values 
of  the  cavitation  index  as  a  function  of  the  flow  coefficient 
at  the  condition  where  the  head  rise  curve  just  began  to 


decrease  due  to  cavitation.   The  details  of  the  impeller 
geometries  are  given  in  Ref.  [6]  and  [8]  and,  for  purposes 
of  application  of  equation  (17),  the  features  indicated  in 
Table  1  were  taken  from  these  references. 


Impel ler    No .  of 


TABLE  1-IMPELLER  INLET  GEOMETRY 
Tan3  at  Vane  Tip 


Vanes   (suction 
side ) 


(pressure 
side) 

(Ave. 

) 

B 

V 

C 

V 

.321 

.306 

.0703 

1.075 

.329 

.314 

.088 

1.095 

.336 

.331 

.105 

1.118 

.46 

.1403 

1.163 

c 

B 
A 
D 


.291 
.299 
.326 


1.04 
1.04 
1.04 
1.04 


It  may  be  noted  that  tanB  in  Table  1  has  two  different  values 
at  the  tip.   This  is  due  to  the  fact  that  the  vanes  were  not  of 
constant  thickness  but  increased  in  width  with  distance  from  the 
leading  edge. 

One  particularly  noteworthy  feature  of  the  experimental  data 
is  the  so-called  velocity  scale  effect.   The  experimental  cavita- 
tion index  in  Figure  3-6  shows  a  decrease  with  impeller  rotative 
speed  with  the  higher  speed  runs  showing  the  better  agreement  with 
the  present  theory.   The  authors  of  Ref. [6]  indicate  that  this 
velocity  scale  effect  has  an  upper  limit,  however,  since  no 
distinct  speed  effect  was  observed  in  a  comparison  of  the  4920 
rpm  runs  with  some  higher  speed  runs  at  6500   rpm.   It  may  be  noted 
that  this  kind  of  dependence  on  speed  has  also  been  found  in  other 
experimental  data  for  submerged  bodies  and,  in  general,  the  velocity 
scale  effect  has  been  found  to  decrease  with  velocity  so  that 
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eventually  K   approaches  1-C    .    .  In  other  words,  the 
o  p  (mm ;  ' 

cavitation  index  approaches  the  value  as  calculated  on  the 

hypothesis  that  cavitation  occurs  when  the  local  pressure 

equals  the  vapor  pressure.   Since  this  classical  hypothesis 

was  the  basis  for  the  present  theory  the  agreement  with  the 

high  speed  runs  is  encouraging. 

The  three  vertical  line  on  Figure  3,4,5  and  6  are  of 

special  interest  when  comparing  the  present  theory  with  the 

experimental  results.   The  one  line  labled  <J>  =tan$  denotes 

the  value  of  <f>   where  the  flow  at  the  upstream  reference 
o 

station  is  such  that  the  relative  flow  angle  with  the  vanes 

would  be  zero.   However,  when  the  fluid  encounters  the  vane 

leading  edges  it  no  longer  has  the  velocity  of  C    but  has 

mo 

velocity  C  .  due  to  the  effect  of  the  shroud.   Thus,  the 
m  1 

condition  of  <j>=tan8(or  d>  =tan8/C  )has  more  significance  than 
1  OS 

<j>  =tan6  since  the  former  indicates  the  condition  where  the 
o 

flow  is  approaching  the  vane  leading  edges  at  zero  incidence. 
However,  as  pointed  out  by  the  authors  of  Ref .  [6J  ,  this  condi- 
tion still  is  not  the  true  condition  of  zero  incidence,  since 
due  to  vane  blockage,  the  relative  flow  just  inside  the  vane 
row  would  not  be  tangent  to  the  vanes.   The  true  zero-incidence 

condition  occurs  when  d)»  =  tan8  or  <i>    =tanB/C  C  .   For  this  condi- 

2  o        s  v 

tion  equation  (17)  reduces  to 


2   2 
K  =C   C 

O    S    V 


which  correpsonds  to  the  minimum  value  of  K  .   This  value  of 

o 

<j>   is  also  indicated  in  the  figures  as  a  vertical  line. 


The  authors  of  Ref.[6j  have  pointed  out  that  their 
experimental  data  shows  that  the  best  cavitation  perfor- 
mance, which  corresponds  to  the  minimum  value  of  K  ,  occurs 

o 

at  the  zero  incidence  condition  of  (J>   =tanB.   Equation  (17) 

corroborates  this  conclusion  to  the  extent  that  K   decreases 

o 

with  decreasing  values  of  A  .  When  $  =  tan6  or  $    =tan6/C  C   the 

u  2  2  °        s  v 

cavitation  index  becomes  K  =C  C   because  the  second  and  third 

terms,  which  arise  on  account  of  the  3  *  /  3 1  term  in  Bernoulli's 

equation,  cancel  each  other.   However,  for  values  of  6    smaller 

o 

than  t  a  n  6  /  C  C  ,  equation   (17)  shows  a  continued  decrease  in 

K   while  the  experimental  data  show  a  minimum  at  d>  =tan8/C  C 
o  o         s  v 

with  a  subsequent  rise  with  decreasing  values  of  <J>  .   This 
continued  decrease  of  K   with  decreasing  values  of  d>   simply 
indicates  that  the  3  § / 3 1  term  occurring  in  Bernoulli's  equation 
gives  rise  to  an  increase  in  pressure  over  that  existing  at 

station  1.   Although  it  is  evident  that  K   could  never  drop 

2 
below  C    (corresponding  to  cavitation  outside  the  vanes  at 
s 

station  1)  experiments  show  a  rapid  increase  in  K   below 

o 

<j>  =tan6/C  C  .   This  rather  extreme  departure  seems  to 

0  s  V 

indicate  that  some  other  mechanism  is  responsible  for  the 

cavitation  breakdown  at  positive  incidence  angles. 

The  difference  between  the  theory  and  experimental  data 

occurring  at  positive  angles  of  incidence,  that  is,  at  values 

of  <i>    =tan8/C  C   is  not  without  an  explanation.   The  one  possible 
o        s  v 

explanation  which  presents  itself  is  that  for  positive  incidence 
the  pressure,  as  calculated  according  to  the  present  theory, 
increases  just  inside  the  vane  row.  (This  is  born  out  by  pressure 
measurements  of  Ref .  \_6j   .)  This  pressure  increase  in  the  direction 
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of  the  flow  (or  adverse  pressure  gradient)  would  naturally 

tend  to  result  in  leading  edge  stall.   However,  for  negative 

incidence  the  present  theory  shows  a  decrease  in  pressure  as 

indicated  by  the  rising  value  of  K   with  <$>     .   This  reduction 

o        o 

in  pressure  inside  the  vane  row  would  provide  a  favorable 
pressure  gradient  and,  consequently,  the  result  would  be  a 
well-behaved,  attached  flow.   Thus,  although  at  the  outset 
there  may  seem  to  be  no  difference  between  positive  and 
negative  angles  of  incidence,  as  would  be  the  case  in  steady 
cascade  flow,  the  8$/3t  term  in  Bernoulli's  equation  gives 
rise  to  an  adverse  pressure  gradient  in  the  one  instance  and 
a  favorable  pressure  gradient  in  the  other.   As  a  consequence, 
the  resulting  flow  is  apparently  altogether  different  for  the 
two  possible  incidence  conditions. 
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SUMMARY  AND  CONCLUSIONS 

A  potential  flow  theory  has  been  developed  to  predict 
cavitation  in  pump  impellers  which  takes  into  account  the 
unsteady  term  occurring  in  Bernoulli's  equation.   The  theory 
has  been  applied  to  mixed-flow  pump  impellers  and  test  results 
show  good  agreement  near  the  zero  incidence  and  overflow 
condition.   As  the  relative  flow  angle  of  attack  increases  the 
theoretical  and  experimental  results  diverge  due  apparently 
to  the  predominance  of  the  leading  edge  cavity  in  causing 
performance  reduction. 

Figure  7  is  an  attempt  to  depict  the  salient  features  of 

cavitation  breakdown  in  pump  impellers.   At  d)  =tanB/C  C   as 

r    r  o        s  v 

indicated  on  the  figure  the  unsteady  term  in  Bernoulli's 

equation  causes  no  pressure  reduct ion, and  equation  (17) 

2   2 
yields  simply  K  =C   C   .   This  point  is  referred  to  as  the 
r  J       o   s   v  r 

point  of  zero  incidence.   For  points  to  the  right  of  this 
point , cor responding  to  values  of  negative  incidence, the  one- 
dimensional  theory, which  disregards  local  effects  at  the 
vane  leading  edge, shows  reasonable  agreement  with  the 
experimental  results.   This  agreement  is  apparently  attri- 
butable to  the  flow  stabilizing  effect  of  the  favorable 
pressure  gradient  caused  by  the  3  $  /  3 1  term  occurring  in 
Bernoull's  equation.   Points  to  the  left  of  the  point 

d>  =tan3/C  C   correspond  to  values  of  positive  incidence, 
o        s  v 

The  sharply  decreasing  K  versus   9   curve  with  increasing 

o  o 

incidence  angle,  a,  indicates  the  existence  of  the  adverse 
pressure  gradients  in  the  inlet  region.   This  adverse  pressure 
gradient  apparently  results  in  leading  edge  stall,  a  phenomenon 
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which  would  tend  to  block  the  flow  passage  and  result  in 
drastic  performance  deterioration. 

Although  the  present  theory  does  not  predict  the  minimum 
value  of  the  cavitation  index  with  acceptable  accuracy  it 
does  shed  some  light  on  a  very  difficult  and  important  problem. 
At  overflow  (or  negative  incidence)  conditions  the  unsteady 
term  in  Bernoulli's  equation  has  been  shown  to  provide  a  very 
important  contribution  to  the  pressure  reduction  in  the  inlet 
region  and  the  theoretical  results  generally  agree  with  the 
experimental  data.   At  underflow  the  time  derivative  term 
contributes  to  an  adverse  pressure  gradient. 
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FIGURE   2.      VELOCITY  POTENTIAL  FIELD 
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FIGURE  3.     CAVITATION   INDEX   VERSUS  FLOW 
COEFFICIENT  —  IMPELLER    A. 
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FIGURE  4.   CAVITATION  INDEX   VERSUS  FLOW 
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